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We analyze the chiral p-wave superconductor in the low temperature region. 
The superconductor has a e x p x + i£ y p y -wave gap in two-dimensional space 



(2D). Near the second superconducting transition point, the system could 

be described by a quasi-lD chiral p-wave model in 2D. The axial anomaly 

occurs in such a model and causes an accumulation of the quasiparticle in an 

inhomogeneous magnetic field. The effect is related to the winding number 
O 



of the gap. 
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Chiral p-wave superfluidity is realized in the superfluid 3 He-A QTJ| . Recently, the possibil- 
ity of chiral p-wave superconductivity has been argued in Sr 2 Ru0 4 . In such superfluids 
or superconductors, the condensation of the Cooper pairs, which have orbital angular mo- 
mentum along the same direction, occurs. Therefore, time-reversal symmetry (T) and also 
parity (P) in two-dimensional space (2D) are violated. We investigate an axial anomaly 
effect in the chiral p-wave superconductor near the second superconducting transition point 
||. It is revealed that the axial anomaly causes P- and T- violating phenomena related to 
the topologically quantized number at sufficiently low temperatures. 

The axial anomaly was originally pointed out in the Dirac QED in 3D @-§[. It is the 
phenomenon in which a symmetry under the phase transformation e* 75 " of the Dirac field in 
the action at the classical level is broken in the quantum theory. Here, a is a constant and 
7 5 is a hermitian matrix which anti-commutes with all of the Dirac matrices 7^, where // is 
the spacetime index. The Adler-Bardeen's theorem guarantees the absence of higher order 
corrections to the divergence of the axial current @||. Therefore, the exact calculation of 
the two-photon decay rate of the neutral ir meson can be carried out. It has been pointed 
out that the same results are obtained by using the path-integral formalism and the relation 
between the axial anomaly and topological quantized numbers has been clarified through 
the Atiyah-Singer index theorem [0 |9[ . 

It has been pointed out that the axial anomaly also plays an important role in the 
quantum Hall effect (QHE) in the 2D massive Dirac QED. In 2D, the mass term of the 
Dirac fermion violates P and T similar to the magnetic field, and the Hall effect may occur. 
It was shown that the existence of the Hall current and its quantization are caused by the 



axial anomaly in ID |T0|- |13|1 . 



The relation between the axial anomaly and QHE in 2D electron gas in the magnetic 
field was also discussed, and the quantized Hall conductance is expressed by the winding 
number of the fermionic propagator in the momentum space ||T2"|-|T6|j . Other applications of 
the axial anomaly to condensed matter physics has been studied in the field of the superfluid 
3 He in 3D USUI, and also in some ID systems H,gO . 



The phenomena caused by the axial anomaly are related to the topologically quantized 
numbers. On the analogy of QHE, it is expected that the axial anomaly also plays an impor- 
tant role in other P- and T-violating 2D systems. In this letter, we show an axial anomaly 
effect in a chiral p-wave superconductor in 2D, which has the spin-triplet e x p x + ie y p y -wave 
symmetry. P and T-violation occur whenever both e x and e y are nonzero. We show that 
near the second superconducting phase transition (2ndSCPT) point the axial anomaly in 
ID causes an accumulation of the mass density of the quasiparticle in an inhomogeneous 
magnetic field. The effect is related to the winding number of the gap; sgn.(e x e y ) fl2~T| , |2"2" 



We use the 2+1-dimensional Euclidian spacetime and the natural unit (h — c — 1) in the 
present paper. 

Let us consider 2ndSCPT in the chiral p-wave state. It has been proposed theoretically 
that such a transition occurs under uniaxial pressure in the x — y plane (the basal plane) 
0. 2ndSCPT caused by an exact in-plane magnetic field has also been argued p3fl , but our 
argument is not applicable to this transition because the anomaly effect which we will show 
requires a magnetic field perpendicular to the basal plane. The orbital dependence of the 
gap near the transition point is shown in Fig. 1. For simplicity, we assume a circular fermi 
surface in the normal state. There are two tiny gap points around p = (±pp,0), where p F 
is the Fermi momentum. Therefore, the quasiparticle excitations along the x-axis become 
dominant in the low temperature region. We can describe such a situation effectively by 
using a quasi-lD model. For the present, we consider the superfluid just for simplicity, and 
will extend our argument to the supercondcuting case by taking into account the Meissner 
effect. We assume a linearized fermion spectrum and a spin-triplet chiral p-wave gap near 
the Fermi surface in the normal state (i.e. p x = ±p? + 5p x , Sp x « p?) written as, 

£r,l(p) = ±v F 5p x , (1) 

A(p) = ia 3 a 2 — (e x p x + ie y p y ) 
Pf 

~ia 3 <j 2 A{±e x + i£ y ^), (2) 

Pf 

where vp is the Fermi velocity. c_r(p) (cl(p)) is the kinetic energy for the right (left) mover. 



e x and e y are parameters and presently they satisfy e x « 1 and e y ~ 1. 
The Lagrangian of our model is written as 



£ = # T 



{^<9 T + n—^ya 3 } <g> 7 T (3) 



A 

-N; F <)p a; l ® 7 X H £ y PyCr 3 ® 7 V - iAe^ <g> 1 

Pf 



*p- 



Here, we use the Bogoliubov-Nambu representation with an isospin a = 1, 2, 



^ tffc(x) ^ 



and \l/ p is its Fourier transform, ^r(x) and V'iXx) are the slowly varying fields for the right 
mover and the left mover with a real spin index s = 1,2, respectively. The matrices 7 r = 
Tlj 7a = 7"2 and 7 y = r 3 are the 2x2 Pauli matrices with isospin indices and Oi (i = 1, 2, 3) 
are the 2x2 Pauli matrices with spin indices. The symbol <jj <S> Jr,x,y represents the direct 
product. \i/ is defined as \i/ = —i^j T . We assume a magnetic field, which is directed to 
the z-axis (the c-axis in the crystal) and has a constant gradient in the ^-direction, i.e., 
B z (y) = (dB z /dy)y, and (dB z /dy) = const. Because we are considering the superfluid for 
the present, the magnetic field couples with the fermion only through the Zeeman term 
nB z ^/<73 eg) 7 r \l/, where /i is the magnetic moment of the fermion. In general, the magnitude 
of the gap A depends on the magnetic field, but we neglect this fact because it has become 
apparent that our result does not depend on A. 



We note that the Lagrangian is equivalent to that of the 2D Dirac QED [24] in a constant 



electric field, except for the existence of 03. It has been shown that the axial anomaly causes 



the quantized Hall effect in the 2D Dirac system [10-13]. In our model, the physical value 
which corresponds to the quantized Hall current j x = e (^bIx^d) i n 2D Dirac theory (\lZrj 
is the Dirac field and < • • ■ > is the expectation value for the ground state) is the mass 
(charge) density, because e (^T^y = e (^r^r + W^l) = e (p e ), where e represents a mass 
of the quasiparticle, and represents an electric charge when we consider a superconductor. 
The reason why e(^ , j x ^) becomes the charge density, not the current, is the remarkable 



difference between \1/d arid \1/. \1/d is a two-component fermion field and each component 
has the same charge. \I/ is also two component (two isospin) but the signs of their charges 
are opposite. 

Let us calculate the expectation value of the mass (charge) density 



<Pe(x)> = Tr 



eo^ 



ifry x + v F 8p x a 3 ® lx - «A£ a 
A hermitian operator d in ID (in the ^/-direction) is defined as 



(4) 



, f „ dB z A A 

d = < lOrO-z + V—T-y^ \ ® It H £yPyl ® 7y, (5) 

where 7 T = -7 y = -r 3 and % = j T = t\. We define 75 as 75 = £7,-% = -7^ = -r 2 , 
and it is anti-commute with 7V and 7 y , therefore, 75 is a hermitian matrix which satisfies 
{75, d} = 0. These findings suggest that if an eigenstate u n of d with a nonzero eigenvalue 
£ n (0 < n) exists (i.e. du n = („«„), 7sM n should be another eigenstate with an eigenvalue 
— £ n . If zero modes of d exist (i.e. du = and ^75^0 = 0), they are divided into two 
groups. One of them is u = (1/2) (1 + 75)^0 with an eigenvalue 75 = +1 and the other is 
u = (1/2)(1 — 75)^0 with an eigenvalue 75 = — 1, since 75 = 1. 

Let us research the eigenmodes of d. The expectation value of d 2 is 

[U n , 4 2 U n ) = \u c \(n+ -) + y(M n ,7 5 M n ), 



dB z 2A 
dy p F 



u c = /j.—^—Ey, (6) 



where u n = u n (y — y c (pT,&3)) is the eigenfunction of the harmonic oscillator with the fre- 
quency uj c . The oscillator is centered at y c (Pr,o'3) = —{dB z /dy)^ 1 {p T /ij)az- Equation (j^) 
indicates that only zero modes which belong to Uq (uq) exist when < u c (uj c < 0). This 
suggests the non-conservation of the vacuum expectation value of the axial charge which is 
defined in the second-quantized formalism as 

(Q 5 ) =(N + -N-), 

N ± = f 'dp y ul (±) u^\ (7) 



while the classical ID theory £id = ^> $$> has the axial symmetry \l/ — > e JQ75 \l/, i.e., the axial 
anomaly occurs. Here, Uq is a second-quantized fermionic field. The anomaly comes from 
the spectral asymmetry of zero modes, as in the discussions presented in Refs. [|I0|- |T3| , |25|| . 
In the free system, the energy spectrum of u Q is po = ± — e y p y for spin down (03 = +1) 
states in Minkowski spacetime, and all of the negative energy states are filled while all of 
the positive energy states are empty and (Q5) = 0. After we turn on the magnetic field 
adiabatically (for the present, we assume < u c ), the energy spectrum of u is lowered 
and (N+) decreases (i.e., empty negative energy states arise on the spectrum of u ); on 
the other hand, the energy spectrum of Uq is lifted and (AL) increases (i.e., filled positive 
energy states arise on the spectrum of u ), and therefore, (Q5) does not conserve. Finally 
(N + ) = and only Uq exists. We see that the same argument is also valid for spin down 
((T3 = —1) states. The nonzero eigenvalues of $ are E n = co c (n + 1/2), since the inner 
product (w n ,75Wn) vanishes whenever $u n 7^ because of the orthogonal relation between 
the eigenfunctions of the hermitian operator. 

Next, we consider the eigenvalue problem of a 2D operator: 



p = fi - AeJ <g> 7^ = fl + A£ x l (g) 75. 



(8) 



Let 



(p n = {a n u n + Pn75U n )e lpx 



(9) 



stand for an eigenfunction. We use a representation for the n-th level such as 

,75 = 



<U o^ 



~ir, 



,U r 



) 



v°y 



(10) 



v 10 ; 



where, 



in 



^\uo c \(n + \) (n=l,2,--.), 
(n = 0). 

Therefore, the eigenvalue equation is written as 



£ n AsJ fa ' 



V 



Ae^ -f„ 



\Pn J 



Cn 



' «„ ' 



\ Ai y 



(11) 



There are two eigenstates for an oscillator in the n(^ 0)-th level written as 



& ±) (p.)=±V£ + A8e 2 



(12) 



a, 






1 
1 



A£„ 



/ 



-Ae x 



/ 



where C± are normalization constants, but for n — 0, there is only one eigenstate 



Co(Px) 



loV 



Wfi 



■A|e x | 



U \ J 



\ Po J 



V 



\ 



/ 



(13) 



-LO c /\UJ c \ 

because the solution should satisfy 75^0 = ~(^>c/\^c\)fo- This condition comes from the 
axial anomaly in the y-direction. 

Finally, we show the accumulation of the mass density from Eq. (HI), which is derived as 



<Pe(x)> = Tr 



E 



-iea 3 



P — iv F 5p x cr 3 Cg> 1. 

dp T f A dSp x 

tr 



2vr J -a 2vr 



e/2dB z y^ r A ddp x 
2vr d V 7^J- k 2?r 



-iea 3 \u n (y-y (p T ,a 3 ) 

CniPx) - W F Sp x (T 3 

1 



-(+)/ 



+ 



^(-)/ 



+ - — / — — tr 



2tt dy J-A 2tt 
sgn(e x e y )efiN 1D (0) 



(n \p x ) - iv F 5p x cr 3 Q \p x ) - iv F 5p x a 3 . 
1 



Co(Px) - iv F 5p x a 3 



dB z 

dy 



(14) 



where the symbol tr represents a trace on the real spin, and we use the normal-orthogonal 
relation / dy\u n \ 2 = 1. A is a momentum cutoff and we assume a relation |A| << A 2 /2m << 
ep. -/Vid(O) = (27TfF)~ 1 is the density of state at the Fermi surface in ID. All of the ti^O 



parts are canceled out because of the co-existence of the eigenvalues (^ and (^~\ Only the 
n = part survives because of the axial anomaly in the y-direction. 

The effect is related to the winding number of the gap in Eq. (|2]) [pl,p2| , 

d 2 p 



16tt 



tr[g ■ (Vg x Vg)] = sgn(ea.e v ), (15) 



' Re[A(p)(-i<7 2 )] ^ 



:(p) 



-Im[A(p) (-icr 2 )] 
y (p 2 /2m) - e F j 

where V = d/dp. This suggests that these effects occur even if e x and/or e y are infinitesi- 
mally small, and that these effects come from the P- and T-violation of the gap. 

The accumulated mass density exists in the bulk region of the superffuid. In the su- 
perconductors, the Meissner effect occurs and the magnetic field cannot penetrate into the 
bulk, therefore the accumulated charge density would exist near the edge |26| of the super- 
conductors and also around the vortex. 

Let us roughly estimate the charge distribution around the vortex under uniaxial pres- 
sure. We make an approximation that the magnetic field is linearly decreased around 
the vortex core and vanishes at r = Al, where r is the distance from the core and Al is 
the London penetration depth; therefore, the gradient of the magnetic field is —H c ij\\^ 
where H& is an upper critical field along the c-axis. We use the parameters appropriate for 
Sr 2 Ru0 4 , H c2 = 0.075 Tesla, Al = 1800 A and vp = 10 6 cm/ s. The induced charge density 
around the vortex core obtained from Eq.fllH]) is of the order of — 10~ 8 e per 1 A 2 , when 
sgn(e a ,e |/ ) = +1. For the charge neutrality, there should be a plus charge n x (1800 [A]) 2 
x 10~ 8 e [At 2 ] = 10~ 1 e at the core. The core charge might extend in the region r < £ for 
the stability. Here £ = ttvf/A is the coherent length. Further (self-consistent) calculation 
is required to obtain more accurate results. 

In summary, we have analyzed the 2D chiral p-wave superconductor near the second 
superconducting transition point |2]|| at sufficiently low temperature. Such a situation 
is described by a quasi-lD model in 2D, which can be mapped onto the 2D Dirac QED. 

8 



Corresponding to the QHE in 2D Dirac theory ]10- |13] , we show that the charge accumulation 



occurs under an inhomogeneous magnetic field as an effect of the axial anomaly in ID. The 



effect is related to the winding number of the gap [£1],^], and could contribute to the vortex 
core charge. 



Recently, the vortex in chiral superconductors has been discussed [ 27| , and such a vortex 
has a fractional charge and a fractional angular momentum. Interesting phenomena related 
to these fractional quantum numbers and the present effects are expected to occur around 
the vortex core. 

An anisotropic chiral p-wave state, which has the d(k) = z(sin k x + i sin k y ) order param- 
eter, is proposed p8| . This state has four tiny gap points, therefore our argument would be 



applicable to quasiparticle excitations in the low temperature region of such a state |29|| . 

The author is grateful to K. Ishikawa, R. Joynt, N. Maeda, K. Shizuya and M. Sigrist 
for useful discussions and encouragement. 
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FIGURES 
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FIG. 1. The Fermi sea in the normal state and the momentum dependence of the gap functions 
for Sr2Ru04 near the second superconducting phase transition point. The shadows show the 
Fermi sea, and the distance between outer lines and inner lines shows the magnitude of the gap 

\A/p F \ 2 (elp 2 x + e 2 y p 2 y ). 
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